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RADIATION FROM SLOTS ON CYLINDRICAL BODIES USING 


GEOMETRICAL THEORY OF DIFFRACTION AND 
CREEPING WAVE THEORY 

3y 

Constantine A. Balanis, Ph.D. 

The Ohio State University, 19&9 
Professor Leon Peters, Jr., Adviser 

A hybrid solution employing wedge' diffraction and creeping wave 
theory is used to compute the principal planes (equatorial and 
elevation) radiation patterns of axial and circumferential slots on 
conducting cylinders of finite and infinite lengths. The slots are 
excited by parallel-plate waveguides operating in the TEM and TE-j_q 
modes. ^ 

For the equatorial plane pattern, the total field in the "lit" 
region is obtained by the superposition of two fields; that is, the 
wedge-diffracted and the creeping wave fields. The wedge-diffracted 
field is obtained by approximating the parallel-plate-cylinder geometry 
with two wedges, each formed by a wall of the waveguide and a tangent 
plane to the cylinder surface at the edge point. The creeping wave 
contribution is obtained by the method commonly employed for the 
computation of scattered fields from curved surfaces. The total field 
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,ih the "shadow" region is obtained solely from the creeping wave 
contribution. For the elevation plane pattern, wedge diffraction 
techniques for the entire pattern are employed. 

j~ The method is checked computationally by comparison with the modal 
solutions for axial and circumferential slots on right circular cylin- 
ders and experimentally for elliptical cylinders since modal solutions 
are not readily available.^ Experimental models are also used for the 
verification of the elevation plane pattern computations since , 
boundary- value solutions are not available for finite length cylinders. 

f Computed results using this technique compare favorably with those 
obtained from existing modal expansion boundary-value solutions and 
experimental resultsjj The main advantages of the present technique are 
that it .can be applied to geometries where modal solutions are not 
possible, in numerical ranges where the convergence properties of modal 
expansions are relatively poor, in parametric design problems since the 
contribution from each field is separated, and in the analysis of 
antennas with finite physical sizes. 
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INTRODUCTION 


The radiation properties of slot antennas on conducting cylinders 
have been studied extensively. s uc h analyses have been restricted 
to geometries where modal expansion field functions are readily available 
and .whose physical length is infinite. Moreover, the solutions obtained 
are poorly convergent for objects more than a wavelength or so in extent. 
However, many practical problems involve slot antennas mounted on finite 
size conducting bodies whose geometrical shape does not conform to a - 
coordinate system where the wave equation is separable. Wedge diffrac- 
tion°~9 and creeping wave theory are used in the analysis presented 
here to compute the radiation patterns of axial and circumferential 
slots on conducting cylinders of finite and infinite lengths. The 
method employed is an extension of the geometrical theory of 
diffraction"^ "^" and creeping wave theory. 

The geometrical theory of diffraction and creeping wave theory are 
extensions of geometrical optics which account for diffraction. They 
introduce diffracted rays in addition to the usual rays of geometrical 
optics. These rays are produced by incident rays which hit edges, 
corners, vertices, or curved surfaces. Some of these diffracted rays 
enter the shadow regions and account for the field. intensity there. 


1 
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The distribution of the incident energy among those diffracted rays is 
described by a diffraction coefficient. The diffraction of the incident 
rays is a local phenomenon, so the diffraction coefficient depends on 
the geometry of the edge or curved surface in the vicinity of the point 
of diffraction and on the polarization of the incident field. 

A field is associated with each diffracted ray and the total field 
at a point is the sum of the fields of all rays at that point. The 
phase of the field on a ray is assumed to be proportional to the optical 
length of the ray from some reference point. Appropriate phase jumps 
must be added as a ray passes through a caustic. The amplitude is 
assumed to vary in accordance with the principle of conservation of 
energy in a narrow tube of rays. The initial value of the field on a 
diffracted ray is determined from the incident field with the aid of an 
appropriate diffraction coefficient. These diffraction coefficients 
are determined from certain canonical problems. 

The canonical problem, which yields the diffraction coefficient 
for an edge, is the diffraction of a linearly polarized plane wave by 
a two-dimensional, infinite wedge. The diffraction coefficient is 
found from Sommerf eld' s 1 ^ asymptotic, high frequency solution to this 

•I Q 

problem. Pauli' 1 ’ has obtained a practical formulation for the diffrac- 
tion coefficient and Hutchins'^ has generalized it. He has obtained a 
series solution which is valid for the exterior and interior regions of 
a wedge. In addition, his solution is accurate for wedges of large 
included angles where the Pauli solution is not accurate. 

Oberhettinger 20 has obtained a similar series in which the leading 

term is identical to the Fresnel integral form for half -plane diffraction. 
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The work of Hutchins tends to bridge the solutions of Pauli and 
Oberhettinger. 

The primary task in applying the creeping wave theory is the 

determination of the diffraction, attenuation, and ray path factors for 

a general body (and it is not yet always practical). Thus it is . 

necessary to evaluate these factors for canonical targets whose exact 

12-l6 21-23 

solutions are available, such as the cylinder and sphere. - 3 ^ 

In order to obtain a more general solution for these factors- it is 

convenient to utilize experimental data to validate an empirical model 

24-26 

for more general targets such as the prolate spheroid. 

The attenuation and diffraction coefficients for a cylinder have 

been obtained by comparing • the asymptotic series of the scattering 

boundary- value solution with that of the diffraction problem 

formulation. 1 -'’" 1 ^ The ray paths along the surface of the cylinder are 

geodesics (great circles). The scattering by a sphere is, however, a 

most important problem. The classical solution of Mie using separation 

of variables and series techniques can be considered as the starting 

21 

point. Senior and Goodrich have obtained an asymptotic form to the 
Mie series solution for the sphere through the application of the 
Watson transformation. Hong has derived attenuation and diffraction 
coefficients with higher order correction terms. The formidable task 
of implementing this solution numerically has yet to be accomplished. 

An approach to the general formulation of the creeping wave paths 
on bodies of revolution has been developed by Kinber. ^ Kinber has 


shown that the wave equation may be expanded in a set of ray coordinates 



in which the solution can he written in a form whose magnitude is 
dependent upon the cross section of the ray tube and whose phase is 
dependent only upon the propagation constant and path length traversed. 

pli 

The empirical approach utilizes a simplified ray path geometry 

to determine approximate attenuation and diffraction coefficients for a 
sphere. This approach is of interest in that it lends itself to 
extension to more general targets. This method is a simplified creeping 
wave analysis. It uses a single nonconvergent (and nondivergent) ray 
path together with approximate diffraction and attenuation coefficients 
to construct an approximate solution for the scattered field due to 
creeping waves propagating in each direction along the ray path. The 
path chosen is the path traversed by the "major" ray (that is, the path 
corresponding to the E-plane of the sphere). This analysis suggests 
that an approximate picture of scattering by a sphere can be constructed 
by neglecting the creeping waves which, have a radial magnetic - field 
(that is, the "minor" creeping waves) and by considering only the 
creeping waves which have a radial electric field (that is, the "major" 
creeping waves) except for the H-plane scattered fields. 

The prolate spheroid represents a body for which an exact closed 
form solution is not available except on the axis of rotation. In 
order to find the ray path geometries for the prolate spheroid for an 
arbitrary angle of incidence, the geodesic corresponding to the point 
of attachment and the tangent direction at that point must be cal- 
culated. The determination of the tangent direction of the creeping, 
wave at the shadow boundary is not an easy task- in general. It is 
therefore suggested that the simplified ray path geometry 2 **"- 2 ^* for the 
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sphere can be applied to the prolate spheroid. In this case, the 
creeping -wave paths are elliptical, thus the total attenuation must be 
expressed as an integral which is dependent upon the radius of curva- 
ture along the path. Also the radius of curvature at the points of 
attachment and reradiation must be computed in order to determine the 

diffraction coefficient. The specific solution for the prolate 

pit. pf^ p^ 

spheroid has been presented by Ryan' * and Peters and Ryan. 

Many other antenna problems have been treated using wedge diffrac- 
tion techniques. Russo, Rudduck, and Peters 2 ® applied the geometrical 
theory of diffraction to calculate the total antenna pattern of a horn 
in the E-plane, including the backlobe region. Theoretical and experi- 
mental patterns are in excellent agreement, thus demonstrating that the 
method for treating diffraction by edges is valid. Obha 2 ^ used the 
geometrical method of diffraction to calculate the radiation pattern 
and gain of a finite width comer reflector antenna. This method 
yields also good results for the computation of backscattering from an 
antenna having conducting plates finite in extent. 

Extensive work of applying the wedge diffraction method to wave- 

6-9 

guides has been carried out by Rudduck. The principal tool employed 
is diffraction by a conducting wedge; the resulting electromagnetic 
field may be treated as a superposition of the geometrical optics field 
and the diffracted field which behaves as a cylindrical wave radiating 
from the edge of the wedges. These techniques may be applied to any 
two-dimensional antenna or scattering body which may be constructed 
from a set of wedges; the radiation or scattering pattern may be obtained 
for any excitation which can be expressed in terms of plane or cylindrical 



6 


waves by superposition of the individual wedge diffractions. Some 
elementary examples include parallel-plate waveguides, walls of finite 
thickness, and polygonal cylinders. 

A basic feature of this technique is that it does not require an 
assumption of the value of the total field or current on some surface, 
as do conventional methods of aperture integration which employ the 
approximation of physical optics. This feature thus allows more 
accurate treatment of problems than ordinarily obtained by approximate 
methods; it also provides knowledge of the fields in terms of the 
incident field, thus resulting in the ability to analytically determine 
admittance and gain of antennas,^ radiation patterns, ^ ^ mutual 

-ZQ 3T 32— 3^ 

coupling, ^ and reflection coefficient. 

The geometrical ray techniques of this method provide conceptual 
simplicity with which solutions may be formulated. This together with 
superposition of wedges allows structural aspects to be taken into 
account, for example, antennas mounted with and without ground planes, 
structures with thin or thick walls, and arbitrary waveguide truncations 
can be treated. In addition, the fields in all space can he determined, 
thus allowing baeklobe regions to he treated. 

Since the introduction of the geometrical theory of diffraction, 
it has been employed successfully in the solution of various types of 
diffraction problems. Previous work in applying the wedge diffraction 
and creeping wave theories, each one individually, to scattering and 
radiation antenna problems has been outlined. However, no one has 
attempted a hybrid solution in using both, wedge diffraction and 
creeping wave theories, oh the same scatterer. The work outlined in 
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the pages to follow, is an attempt to formulate a hybrid solution 
making use of both wedge diffraction and creeping wave theories on the 
same' antenna. The antenna of interest is a slotted, axially and 
circumferentially, conducting cylinder. The slots are excited by 
parallel-plate waveguides operating in the TEM and TE 10 modes. 

The essential feature of this approach is that it is applicable to 
all .types of complicated problems, some of ^which are impossible to solve 
rigorously. The approach is to resolve a complicated problem into 
simpler ones, each of which will have relative simple rigorous solu- 
tions. The. basic building blocks will be rigorous solutions to such 
canonical diffraction problems as the wedge, cylinder, sphere, etc., 
which will he put together to solve more complicated problems. 

The work outlined in this hybrid solution will .be for general 
conducting cylinders of arbitrary convex cross section. However, com- 
putations will be restricted to circular cylinders for which boundary- 
value solutions exist for comparison and elliptical cylinders for which 
models and experimental data are available for' comparison. The radiation 
patterns of interest for each case are the principal plane patterns, 
equatorial plane (0 = 90 °) and elevation plane (0 = 0°) in the usual 
spherical coordinate system. 

The equatorial plane pattern will he obtained by the superposition 
of diffracted fields from a set of wedges and creeping wave fields 
propagating around the surface of the cylinder. The set of wedges 
approximate the parallel -plate-cylinder geometry in the vicinity of the 
discontinuity formed by the waveguide and the cylinder. Each wedge is 
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formed by a. wall of the parallel-plate waveguide and the tangent plane 
to the cylinder surface at the edge point. 

For the equatorial plane pattern, two distinct regions can be 
identified: the "lit" and "shadow" regions. The "shadow" region 

encompasses the space occupied by the set of wedges and cylinder surface 
and the "lit" region the remaining space. The total field in the "lit" 
region is obtained by the superposition of two fields; that is, the 
wedge -diffracted and creeping wave fields. The total field in the 
"shadow" region is obtained solely from the creeping wave contribution. 
For the elevation plane pattern, wedge diffraction techniques for the 
entire pattern are employed. 

The most interesting part of this hybrid solution is the coupling 
mechanism between the wedge-diffracted and creeping wave fields. 

Creeping waves are launched by the tip of the wedges in a direction 
tangent to the cylinder surface at the edge point. Since the field 
must be continuous along the line separating the "lit" and "shadow" 
regions, the wedge-diffracted fields along the surface of each 
wedge will be used as the initial value of the creeping waves to 
preserve the continuity of the fields. The creeping waves continually 
shed energy in a tangential direction as they propagate around the 
cylinder surface until they reach the opposite w all of the slot and 
again illuminate the wedges; this is a second order field contribution 
and it is neglected. The loss of energy due to reradiation is accounted 
by the use of an attenuation factor. 

The method is checked computationally by comparison with the modal 
solutions for axial and circumferential slots on right circular cylinders 
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and experiment al l y for elliptical cylinders since modal solutions are 
not readily available. Experimental models are used for the verifica- 
tion of the elevation plane pattern computations since boundary-value 
solutions are not available for finite length cylinders. 

Computed results using this technique compare favorably .with those 
obtained from existing modal expansion boundary-value solutions and 
experimental results. The main advantages of the present technique are 
that it can be applied to geometries where modal solutions are not 
possible, in numerical ranges where the convergence properties of modal 
expansions are relatively poor, in parametric design problems since the 
contribution from each field is separated, and in the analysis of 
antennas with finite physical sizes. 



CHAPTER I 


WEDGE DIFFRACTION 

Wedge diffraction techniques are applied for the analysis of the 
radiation pattern of the parallel-plate waveguide shown in Figure 1. 
This geometry has been treated using wedge diffraction concepts by 
Rudduck and his co-workers. Their .solution is incorporated in the 
slotted cylinder radiation pattern analysis and. will be outlined below. 
In this approach a single-diffracted wave emanates from each wedge that 
is illuminated by an incident plane wave. The single-diffracted waves 
produced by one wedge illuminate the other producing double-diffracted 
fields. This process continues to higher orders of diffraction. 

Single Diffraction 

The principal method employed in the analysis of a parallel-plate 
waveguide is diffraction by a conducting wedge. The diffraction of a 

1 y l8 

plane wave by a wedge was solved by Sommerfeld. Pauli obtained a 

practical formulation of the solution for a finite angle conducting 

19 / \ 

wedge which was later improved by Hutchins (see Appendix). The 
total electromagnetic field from the wedge may be treated as the super- 
position of geometrical optics and of diffracted fields which behave as 
cylindrical waves radiating from the edge of the wedges. 
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Figure 1. Geometry of a parallel-plate waveguide aperture. 






The diffraction of a plane wave by a wedge is shown in Figure 2. 
The solution to the plane wave diffraction problem may be expressed in 
terns of scalar functions that represent the normal and parallel 
polarization components of the electromagnetic field to the plane of 
study in Figure 2. The total field is defined as 

E = E g + Ep (D 

where Eq is the geometrical optics field and Ep is the diffracted 
field. The diffracted field is given by 

% = V B (r,i|r - ^o> n ) * v B^ r ^ + 

where the parameters r, ty, n are defined in Figure 2 and 

Vg(r,i]r ± n) is the diffraction function as defined in the Appendix. 
The plus (+) sign applies for the polarization of the electric field 
normal to the edge 



wedge 


= 0 


(3) 


where n is the unit normal to the edge of the wedge and the minus (-) 
sign applies for the polarization parallel to the edge 


Iwedge) ® ^ 

The geometrical optics field is defined in three regions, as shown in 
Figure 2. For plane wave incidence, the geometrical optics field is 
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Reflected Geometric 

Incident Geometric Field ± Field Regions 


e G = e 


ykrcosOMO 


rt - i 0 <i!f<it + t 0 (5) 


jkrcosCt-^n) , jkrcos(t+t 0 ) 
Ep = e u ± e 


0 < t < jr - (6) 


% = ° 


it + t 0 < ^ (7) 


The time dependence is used throughout this analysis. 

The diffracted wave Ep may he represented as a cylindrical wave 
radiating from the edge (see Appendix). In fact, at large distances 
from the edge and in regions removed from shadow boundaries Ep has 
the radial dependence e~^^ oc / Vr. Because of this cylindrical nature, 
subsequent diffractions of a diffracted wave may be treated as the 
diffraction of a cylindrical wave by a wedge. 

Diffraction By a Pair of Wedges 

The process of diffraction by a pair of wedges (parallel-plate 
waveguide) w ill be outlined below. Two polarizations which must be 
distinguished are those corresponding to the TEM and TE^q waveguide 
modes. The two cases of interest are shown in Figure 3. In the TEM 
mode, the incident plane wave is parallel to the axis of the guide 
having a polarization perpendicular to the edge of each wedge. The 
TEpo waveguide mode may be represented by two TEM waves which reflect 
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obliquely back and forth between the waveguide walls. ^ Thus, in the 
TEio mode the incident plane waves have a polarization par all el to the 
edge of the wedges. 

1. TEM Mode 

In the TEM mode an incident plane wave propagates parallel to the 
axis of the guide walls as shown in Figure 3(a). The far-zone singly 
diffracted fields from wedges 1 and 2 are given by 



The superscripts denote the order of diffraction. Expressing the phase 
of the incident waves with respect to edge 1 


E^ = E = 1 


( 10 ) 


E 


2 


Ee -3kwcot0 g _ e -jkwcot0g 


(11) 


Equations 8 and 9 reduce to 
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E Dl^( r l , ^l) “ 


-3 [kr*i +'*A] 

\j2n'kr i 




( 12 ) 


E^2^V ^ 2 } = 


where K^^) and R^ ; (0 2 ) are the singly diffracted rays from 
wedges 1 and 2 given by ' 


3 -j(hr 2 +*/i0 
\j 2jrkr2 




( 13 ) 


>( 1 ), 


- 




- jkwcot0 


g i 
n 2 


frt 


COS' — — cos 
n 2 ' 



( 14 ) 


( 15 ) 


Applying the far-field approximations 

h ^ 

r « r + ^ cos(0 + 0 ) 

1 o 2 o r g' 

) 

r 2 ~ r o - | COS (0 O + 0g) 
r l ~ r 2 “ r Q 


For phase terms (16) 


For amplitude terms (17) 
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where 


h “ *2 ” 

lead to 


(18) 




-j {k(r o +(h/2) co s(0o + 0g3 + «A} 

\/2jrkr 0 


41 ’ 


( 19 ) 




-att? o +(h/ 2 )eos (0 o +0g )] + *A } 

e 




( 20 ) 


The singly diffracted rays may again he diffracted producing doubly 
diffracted rays and so on to higher orders of diffraction. The singly 
diffracted ray illuminates edge 2, giving rise to the doubly 

diffracted ray and R^ causes in a similar manner. 

Also, some of the singly diffracted rays from edge 1 will be reflected 
from wedge 2 and appear to radiate from the image of edge 1 giving 
rise to 

. 4 & 0 O > - (21) 

The doubly diffracted ray from edge 2 is given by 
R d1 } (0o) = E D1G&bA^ " 0 O ~ 0g^ 2 ) + V bA'* ~ 0 O + n 2 0 ^ 



19 


where 

4lG = Sdi («) 

and from edge 1 "by 

SkEX* - + *g'”l> + V 11 ’ 2 * + K ■ *V n l>] (24) 

where 

4ld - 4a* (* - S* g > (25) 

There is an. additional doubly diffracted ray by edge 1 from the image 
source which is given by 

E D1R^o) = + n lj + V b( 2w '^ + ^ n l) ( 2 ^) 

where 

- 4i{~ f ) W 

The total diffracted rays from edges 1 and 2 , using single-double 
diffractions, are given by 

SraW " + 4i ( ^> + 4irWo> ( 28 ) 

e D2^o) -41’^d) + 4a > (* < o> (29) 
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and the reflected ray from wedge 2 by 

%FL^o) = R Dl("0o) ^ ° 

The total diffracted field from the aperture may be expressed as 
the superposition of the total diffracted rays from edges 1 and 2 plus 
the total rays from the image source to yield 




-3 (kjr o +(h/2)cos(0 o +0 g )j +*A } 

e 

\J 2irkr 0 



R D2^o' 1 


JKhcos(0 o +0 ) 




■j2kwsin0 o ^ 


(3D 


Each tern in Equation 31 contributes to the radiation pattern only in 
certain regions as follows: . 

•s 

Boxtfo) -fig < i a < « - WA1 (32) 


B D2 (0 o ) k + WAS < 0 O < it - fi g (33) 

Bpif-fo) \ < *o<l (34) 

There are subsequent diffractions which result in triple and higher 
order diffractions from edges 1 and 2. The total higher order diffrac- 
tions (that is, second and higher order) can be put in a closed form. ^ 
The total illumination of edge 2 from edge 1 can be expressed as 
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r ig = R i(-^ g ) (35) 

where R]_(0 O ) is the total diffracted ray from edge 1. Consequently, 
the total higher order diffractions from edge. 2 are 'given hy 

42 )( <*o> ' S l G [V h ”‘ - t + V h ’* ' K + ^” 2 0 ' (36) 

and the total ray by 

E #o> “ + E D2 ) (i'> (37) 

The total higher order illumination of edge 1 is given by 

E 2G = R 2 ( * ~ 0g } ’ (38) 

and 

E 1P - E i(- f ) (3?) 


Thus the total higher order diffractions from edge 1 are given by 


Sdi’gU - + 0 g *”i> + v B ( h . 2 ” + - (^“lO 

+ B lp jv B (2»-,f + j, 0 , n x ) + Y B (2w,^ + 
and the total diffractions by 

E l(*U - KdI^o) + 4?(0 


(HO) 


0*1) 
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The total diffracted field from the aperture can he expressed as 
the superposition of the total diffracted rays from edges 1 and 2 plus 
the total reflected ray as given hy Equation 31 where %) 1 (0 O )^ 

%2^ 0 )> ^ r di(" 0 o ^ ^ re P laced ^y %(0 O )> ^(00^ E i(-0 O )^ 

respectively. 

The total diffracted waves from edges 1 and 2 are given in terns 
of the unknown illuminating rays Rpg, ^2G’ R 1P* These rays can 

he determined hy the solution of three simultaneous linear equations 
formed hy expressing each unknown ray in terms of Equations 37 or 4l. 


R 1G 





Rip = R- 


(i) 

1 D1P 


+ RonV' 


2G V 2G 


(- f ) %>%{- f) 


m 


%} “ E D2G + R 1G V 1G^ ~ ^ 

where the quantities V__, V__, and V are the unit-wave 

Xtr IP PLr 

diffractions used in Equations 36 and 40. 

2. TEpo Mode 

The diffraction at the aperture of the parallel-plate waveguide 
for the TEp 0 mode may he treated in a similar manner as the TEM mode. 
The TEpQ mode may he represented hy two plane waves reflecting 
obliquely hack and forth between the waveguide walls at an angle 
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= sln ' 1 (^) (45) 
as shown in Figure 3(h). The wave has a polarization parallel to the 
edges of the wedges which- form the waveguide walls. For this mode 
two cases must he distinguished, as determined by 


Case I: A > 0 (46) 

o r g 

Case II: A Q < 0 g (47) 


In Case I edge 2 is not illuminated "by the incident plane wave and no 
singly diffracted ray emanates. For Case II both edges are illuminated 
hy the incident wave; hence, singly diffracted rays emanate from both 
edges. 

Since the polarization is parallel to the edges -of the guide for 
the TEpo mode, the negative (-) sign in the. diffraction formula is 
chosen. The singly diffracted ray from edge 1 is obtained from 
Equation 2 as 


Kel’^o) = 



COS 


n l 


- cos 


fa + K k o 


n l 


It fa + 0o + A o 

cos cosl 

n l \ n l 


(48) 


For Case II (A^ < 0_), the singly diffracted ray from edge 2 is given 
U — & 
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iui.(0 ) = - — sin — e 
V2 sr o' n 2 n 2 


- jkw jsinA o +cot0gCos 



(49) 


where the exponential factor and minus sign represent the phase of the 
incident plane wave at edge 1. The reflected rays are given by 

*sn,V o> = h < t < f (50) 

where the preceding minus sign results from the reflection. 

Multiple diffractions occur in the same manner as for the TEM mode 
but with the minus sign chosen in Equations 22, 24, 26, 36, and 40. 

Thus the total higher order diffracted waves from edges 1 and 2 for 
the TE 10 mode are given by 



and 


- K - - V 1 ** - K - k'S 


(52) 
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The minus sign preceding R-jj results in the same manner as for the 
reflected rays of Equation 5 0 • The total wave from each edge is 
obtained by using TE-^q mode rays in the same equations valid for the 
TEM mode, that is. Equations . 57 and 4-1-44-. The unknown illuminating 
rays are determined in the same manner as for the TEM mode by using 
the formulations for R^, R^* R^, and given above for 

the TEpo mode. 



CHAPTER II 


CREEPING WAVE DIFFRACTION CONCEPT 


When a wave is incident upon an opaque object which is large 
compared to the wavelength* a shadow is formed. However* some radia- 
tion penetrates into the shadow region due to diffracted rays d ^ as 
shown in Figure 4. These rays are produced by incident rays which are 
tangent to the surface of the body. Each tangent ray splits at the 
point of tangency with one part continuing along the path of the 
incident ray and the other traveling along a geodesic on the surface 
of the body. At each following point* it splits again with one part 
traveling along the geodesic and the other reradiating along a tangent 
to the geodesic. From a single incident ray, infinitely many diffracted 
rays are produced* one of which is reradiated at each point of the 
geodesic. These waves traveling around the opaque body have been 

12 

designated as creeping waves introduced first by Franz and Depperman 
for the interpretation of scalar diffraction by circular cylinders and 


spheres. 

The scattered field caused by the creeping wave mechanism for a 
plane wave incident on a cylinder is given by 


Ep(s*0) = E.(Q) 


5 -Jk(t+s) 

VS 


I 


m 




mh 


(P)e 


- f Omh(p) ds 


(53) 
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■which for a right circular cylinder reduces to 


Ep(s,0) = E,(Q)- 


-jk(t+s) V -1 2 ~ a li]h t 


'/s' 


l 


m 


D mh e 


( 54 ) 


since the radius of curvature is constant. 


D mh = D tri -h (Q) D mh ( F ) = Diffraction Coefficients 


(55) 


< x mh = Attenuation Constant 


~(kp) 

2p 


1/5 


2/5 


^(4m + N) 


,3 *, 


l<5 


(5 6 ) 


p is the radius of curvature of the body’ 
k is the propagation factor 
t is the path length along the body 

s is the distance from point of detachment to the observation point 

P is the point of detachment of the creeping wave 

Q is the point of attachment of the creeping, wave 

For a hard surface (E-field normal to the surface) d Q = 1.083, 

d-j_ = 0. 555, N = 1, and for a soft surface (E-field parallel to the 
surface) d Q = 0.645, dj_ = 0.4-90, E = 5- The diffraction coefficients 
and attenuation constants are obtained by comparing the asymptotic 
series expansion for large bp of the canonical boundary-value problem 
with Equation 55. !5>l4 
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Slotted Cylinder Creeping Wave Mechanism 

Rays diffracted by the edges of a slotted cylinder shown in 
Figure 5 which are tangent to the surface of the cylinder hehaye in a 
similar manner. At each apex point additional diffracted rays are 
introduced each of which will travel along the surface of the cylinder 
and will reradiate' tangentially. The creeping wave field for a slotted 
conducting cylinder can be. expressed as 


"jk(tp+Sp) r-1 - I 

E c (s,0) = %(Qp)- — ) ^mh ( %. )®mh (^*1 ) e 

& m 

-Jk(t2+S2) V -1 

+ Ei (Q2)^ — ) D n j 1 (Q 2 )D nh (P2)e 

V®2 • 


%ih(p) ds 


' Onh(p)ds 

° (57) 


n 


Considering only the lowest-order creeping wave mode. Equation 57- 
reduces to 


E c (s,0) - E i fQ 1 )D oh (0. 1 )D oh (P 1 )- 


E i (%,)D oll (q 2 )Doii(P 2 )^ 


_ f\ oh ( p)ds 
e Jo 

V' 5 ! 

pt 2 

-jk(t2+S2) _“J 0 a oh(p) ds 


(58) 


which for a circular cylinder simplifies to 


E C (^0) = E i (Q 1 )D^ h (Q 1 )^ 


Ei(Qe)D| h (Q2)- 


2 , -ak(ti+r) -aohtq 


r jk(t 2 tr) t 


'oh 2 


s/t 


(59) 


since the local radius of curvature is the same at Pp, P 2 , Qp, and Q 2 . 
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To assure continuity of the fields across the boundary between the "lit" 
and "shadow" regions, the far-field edge diffracted fields are to be 
equated to the creeping wave fields at the shadow boundary. The 
angular variation of the wedge-diffracted fields along the shadow 
boundaries is given by 




for the geometry of Figure 6. In other words, the angular variation of 
the wedge- diffracted fields along the shadow boundary will be used as 
the initial value of the angular variations of the creeping waves. 

Thus, and Eg will serve as the coupling mechanisms between the 

'wedge-diffracted and creeping waves which assure field continuity along 
the shadow boundaries. The creeping wave field for a circular cylinder 
can then be expressed as 


En(2\»0) — — 


jkCt-wtr) -jk(t 2 +r) 

e^l + R B «■■*■■■- e~ a< 
vr \/v 


’oh.^2 


(62) 


where 


r A ®i^l^oh^l^ 


e -pkr 




( 65 ) 
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% - - (| + e) 


For a noncircular cylinder the diffraction coefficient D oh is a 
function of the radius of curvature. 'Therefore the creeping wave field 
given by Equation 62 must be multiplied by D 0 ] a (p)/D 0 ] 1 (Q,) giving 

P t 1 , . 


E c (s,0) = E a 


D oh( p l) e -jk(t 1 +s 1 ) ^ f Q a ° h P 
D oh^l) 

Pt 2 . 

-«ik(tp+sp) J a oh^ p ^ d " 


B D oh ^) \/®2 


The wedge approximation of the parallel-plate-cylinder geometry 
does alter the physical boundaries of the antenna structure and the 
diffraction mechanism of the fields especially in the penumbra region. 
For the wedge approximation, the fields in the 0 = ±(jt/2 +£) directions 
travel along the boundaries of the wedges. However, in the actual 
antenna structure some of the energy is trapped on the surface of the 
cylinder and travels along the curved surface reradiating in a tangen- 
tial direction. Thus the wedge-diffracted and creeping wave fields in 
the transition region (penumbra) will not satisfy the field boundary 
conditions. However, computations carried out using such a model com- 
pare favorably with existing boundary-value solutions and experimental 
data. Therefore, such a model would be a valid approximation of the 
structure as far as the diffraction mechanism of the fields is concerned 
and Equations 60 and 6l would serve as the coupling mechanisms between 
wedge-diffracted and creeping wave fields. 



CHAPTER III 


CIRCULAR CYLINDER RADIATION 

The wedge diffraction and creeping wave techniques previously 
outlined will first he applied for the calculation of the equatorial 
radiation patterns of axial and circumferential slots on circular 
conducting cylinders of infinite length. Boundary-value solutions 
for slots on circular cylinders with common feeds exist and will he 
used for comparison. Once the method is verified, it will he employed 
for pattern calculations of slots whose field distribution is such that 
modal solutions do not exist and to bodies whose- geometric shape does 
not conform to -a coordinate system where the wave equation is separable. 
In addition, the elevation plane pattern for finite length cylinders 
will be analyzed. It should be pointed out that no boundary-value 
solutions for finite length cylinders exist, and experimental results 
will be used for comparison. The diffraction contributions from the 
edges of the cylinder and the effect of' the finite width aperture to 
the overall pattern will be observed. Pattern calculations for more 
complex geometries such as an elliptical cylinder will be carried out 
in the following chapter. 

Equatorial Plane Pattern of Axial Slots Operating in the TEM Mode 

The two -dimensional diffraction geometry for an axial slot mounted 
on a circular cylinder and operating in the TEM mode is shown in 
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Figure 7 • The junctions formed by the walls of the parallel -plate 
waveguide and the planes tangent to the surface of the cylinder at the 
edge points are represented by a pair of infinite wedges of finite 
included angle WA = (2 - n ) n. Two regions, "lit" and "shadow," are 
formed by the imaginary sides of the finite wedges as shown in Figure 7* 
For 0g = 90° of Figure 3(a), K D1 (-0 O ) = 0, %lr( 0 o ) = 0, 
h = w, = ng = n and Equation 31 for the wedge -diffracted field 
reduces to 


E D( r oA) = 


- j { k (r 0 - (w/ 2 )sin0]' + it fk } 
\j 2irkr 0 


%(0 O > + E 2 (0o)e' 3kWSin00 ] 


( 66 ) 


A shift of coordinates from the aperture to the center of the cylinder 
is convenient, so that common coordinates will be used for wedge- 
diffracted and creeping wave fields. Assuming the fan -zone approximations 


r ~ r - a cos 0 cos 3 


r o * r 


00-0 


For phase terms ( 67 ) 


For amplitude terms (68) 


( 69 ) 


the total wedge-diffracted field is- expressed as 


E D (r,0) 


1 

\/2rtk 





3 -jkwsin0^ 


(TO) 
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The creeping wave contribution as expressed by Equation 62 is 
f given by 


E n (r,0) = R 


e -jkr -t^a^+Jk) e"^ kr _t 2^ a ’oh + ^ k ) 

‘A — e + e B -ft- e 


(71) 


The wedge-diffracted and creeping wave fields, suppressing the 
e- 31a 7 \Tr factor, are given by 


% (0) . -Jg e 3 { E &/2)sin0-teeos0cos P ] -«A) 7^(0) + 

(72) 

and 

E O (0) = + (75) 

In general, the total field is equal to the sum of the two fields 

E t (0) = E d ( 0) + E c (0) . (74) 

However, wedge-diffracted and creeping wave fields do not exist in all 
regions. To find the total field in each region, the appropriate wedge- 
diffracted and creeping wave fields must be considered, as they are 
tabulated in Table 1. 

Radiation patterns computed using the boundary-value solution given 
by Wait ^ and the fields from Table 1 in their respective regions for the 


diffracted-field solution are shown on Plates I, II. It should be noted 
that the second and higher order diffractions are approximated by 



TABLE 1. WEDGE-DIFFRACTED AND CREEPING WAVE FIELDS FOR THE EQUATORIAL PLANE PATTERN 
OF A CIRCULAR CYLINDER IN THE DIFFERENT REGIONS 
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wedge diffractions resulting from uniform cylindrical waves. Diffracted 
waves will resemble cylindrical waves if they are observed sufficiently 
far away from any shadow boundary that results from their source. The 
double-diffracted waves can be adequately treated as those of uniform 
cylindrical waves for guide widths down to about A/5* However, triple 
and higher order diffractions are not very accurate -when 0 g = 90° 
since they are viewed at the shadow boundary formed by their source^ 
and the solution is more accurate for wider slots. 

The solutions which consider only first and second order diffrac- 
tions have a discontinuity in the 0 = ±90° directions because higher 
order diffractions are neglected. The discontinuity and variation from 
the boundary- value solution in the 0 = ±90° is greater for smaller 
guide widths. The solution which takes into account higher order 
diffractions eliminates the discontinuity, but it is ,not very accurate 
in the geometrical shadow boundary region (penumbra) since the assump- 
tion of uniform cylindrical waves is not satisfied for triple and 
higher order diffractions. 

This analysis is more accurate for large 'ka cylinders since the 
decaying exponents of the surface waves are derived from asymptotic 
series for large ka. Using this approach, the computations for large 
ka cylinders do not require any additional details from those of small 
ka. In contrast, the boundary-value solution computations for large ka 
cylinders are not very convenient since Bessel functions of large order 
and argument are not readily available. Also, the convergence of the 
radiation pattern function becomes poor for large ka, and higher order 
terms must he included (about 2ka terms). From an economic point of 
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view, large ka value computations using the boundary-value solution 
require considerably more computer time compared to the corresponding 
diffracted field solution (about a factor of 50) since higher order 
terms must be included for convergence. 

As the values of ka increase, the field in the "shadow" region 
-(jt/2 + p) > 0 > (jt/2 + (3) decreases and the ripples start to appear 
at larger angles. It is noted that the period of the ripples is 
approximately (l8o/ka) degrees. A field plot (voltage versus angle) for 
the ka = 12, w/A = 0.2 cylinder is shown in Figure 8. The variations 
between the boundary-value and wedge-diffracted solutions and the field 
discontinuity along the 0 = ±90° directions are more noticeable on a 
linear scale. The agreement would be significantly improved for larger 
cylinders with wider plate separations as shown for the ka = 36, 
w/A =0.4 cylinder in Figure 9- 

Up to now, the parallel-plate-cylinder geometry has been approxi- 
mated by a set of wedges of finite included angle (n ^ 2) each formed 
by a wall of the parallel-plate and a tangent plane to the cylinder 
surface at' the edge point. For the far field, it would be the presence 
of the conducting wall for the wedge geometry of Figure 3 that would 
force the 'field to be zero at 0 = it - WA. However, there is no such 

conducting wall in the far field of the geometry of Figure 6 and it is 
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only the edge that needs to be considered. Sommerfeld 1 points out that 
patterns on precise diffraction photographs exhibit almost no dependence 
on the material and shape of the diffraction edge, and a glass surface 
with radius of curvature of several meters yields essentially the same 
diffraction fringes as the edge of a razor. It is then suggested that 



Diffracted-field (multiple diffractions! 
Diffracted-field (single-double diffractions) 
Boundary-value 



270 ° 


180 ° 


Figure 8. Field plot of axial infinite slot on circular conducting 
cylinder (ka = 12, w/A = 0.2) using finite wedges (TEM mode). 




Diffracted-field (multiple diffractions) 
Diffracted-field (single-double diffractions) 
Boundary -value 



1270 ° 


Figure 9* Field plot of axial infinite slot on circular conducting 
cylinder (ka = 36, w/l = 0.4) using finite wedges (TEM mode). 
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another possible approximation of the parallel-plate-cylinder geometry 
will be a half -plane (wedge with n = 2) instead of a finite wedge 
(n f 2). This approximation becomes necessary for the TEpo mode because 
the wedge-diffracted field along the boundary separating the "lit" and 
"shadow" regions for the finite wedge approximation will be zero. Thus, 
there will be no creeping waves traveling around the cylinder surface 
for the TEpo mode and no fields in the "shadow" region. Since this is 
not true, the half -plane approximation becomes necessary. 

However, to verify the validity of the approximations, finite 
wedge (n ^ 2) and half -plane (n = 2), both models will be used for the 
computation of patterns for TEM and TEpo mode slots and compared with 
boundary- value solutions and experimental results. 

It was found, by comparison, that another approximation for the 
TEM mode slot on a circular cylinder was to replace the edge 1 geometry 
by a half -plane (n = 2) for the first order diffraction and by a finite 
wedge (n ^ 2) for second and higher order diffractions for' the 0°-l80° 
pattern measured in the counterclockwise direction. The edge 2 geometry 
was replaced by a finite wedge for all orders of diffraction. For the 
l8o°-36o° pattern, the approximations of edge 1 geometry are valid for 
edge 2 and vice-versa. 

The computed results using the above approximations are shown on 
Plates III, IV where they are compared with the boundary-value solution. 
It is seen that the above approximations result in better accuracies for 
smaller guide widths around the penumbra region than the finite wedge 
approximation used previously. However, the finite wedge approximation 
gives better results in the "lit" region as it should since the fields 



Plate III. Radiation patterns of axial infinite slot on smaller 
circular conducting cylinder using half -plane and finite 
wedges (HM mode). 



Plate IV. Radiation patterns of axial infinite slot on larger circular 
conducting cylinder using half -plane and finite wedges (TEMmode). 
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in the "lit" region are more strongly dependent on the geometry of the 
edge of the wedge, while the far fields in the penumbra region would be 
dependent on the fields on the walls of the wedge. It should then be 
concluded that both approximations are valid and a combination of the 
wedge approximations dependent on the location of the observation point 
would yield optimum results. As a first order of approximation, either 
one of the suggested models may be used. A field plot for the ka = 12, 
w/A = 0.2 cylinder is shown in Figure 10. The accuracy of the above 
approximation of the penumbra region is better seen in this plot when 
compared with the one in Figure 8. 

As was pointed out earlier, many models for the parallel-plate - 
cylinder geometry would give good results as a first order of approxi- 
mation. Another suggested model for the parallel -plate-cylinder 
geometry would be a thin-walled waveguide (n = 2 for all orders of 
diffractions). Computed results using such a model are shown in 
Figure 11. A reasonable agreement between the boundary-value and 
diffracted solutions is indicated. 

Sometimes it may be desirable to have the maximum radiation 
oriented at a given angle other than the zero degree direction 
described previously. One way of accomplishing this would be to have 
the feed of the waveguide mounted in a slanted position. 

TEM mode propagation along the parallel-plate is assumed but the 
aperture field distribution will not be uniform in this case. Propaga- 
tion of this mode in a guide mounted on a cylinder in a slanted position 
has no boundary-value solution. However, the technique of diffracted 
fields described previously can be used to obtain the radiation pattern. 



hs 


• Diffracted-field (multiple diffractions) 
Diffracted-field (single-double diffractions) 
. Boundary-value 



180 ' 


Figure 10. Field plot of axial infinite slot on circular conducting 
cylinder using half -plane and finite wedges (TEM mode). 
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Diffracted-field (multiple diffractions ) 
Diffracted-field (single-double diffractions) 
Boundary-value 



270 ° 


ka = 12.0 
1 ^ 0.2 


180 ° 


Figure *11. Radiation pattern of axial infinite slot on circular 
conducting cylinder using half -planes (TEM mode). 


The two-dimensional geometry of this mode is shown in Figure 12. 


The modified equation for the diffracted field is given by 


E l c d (k[acos0cosp-(h/2)cos(0+€+0 a )] - jiA) 

\/2rtk 

-^Kl(0) 


. „ (r i, jkhcos(0+€+0 ) . -O2kwsin(0+e} > 

+ Rg(0)e g + 

j 

1 

(75) 

where 



% 2 


(76) 

h "“V 

sin- 0 g 


(77) 


and the regions where the appropriate wedge-diffracted fields are 
applicable are given by 


Rq(0) 

l+ p >0>_! 

(78) 

R 2 (0) 

§ > 0 > - (f + ?) 

(79) 

^(-0) 

— - e > 0 > (0 - e) 

2 . g 

(80) 


The creeping wave field is the same as that given by Equation 73 
applicable in the regions outlined in Table 1. 
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The computed radiation patterns for this mode of operation are 
shown on Plate V with a field plot for one cylinder in Figure 13. It 
is observed that the maximum radiation for the cases investigated is 
oriented at an angle 0 = -e/2. This agrees "with the physical reasoning 
by observing the reflection mechanism of the rays from edge 1 by the 
side of wedge 2. It can then be concluded that by slanting the feed 
waveguide at a given angle, it will orient the maximum radiation in 
other than the 0 = 0 ° direction. The discontinuities observed in the 
solution which considers only single-double diffractions occur because 
third and higher order diffractions which are neglected become 
significant. 

Equatorial Plane Pattern of Circumferential Slots Operating in the 
TE 10 Mode 

The TE i0 mode propagation in a parallel-plate guide can be repre- 
sented by two plane TEM waves reflecting obliquely back and forth 
between the waveguide walls with the electric field parallel to the 
edge of the walls. The approximation of the parallel-plate- cylinder 
geometry by a set of wedges each formed by a waveguide wall and a 
tangent plane to the cylinder surface at the edge point will not satisfy 
the electric field boundary conditions. In order to overcome this 
obstacle, the half -plane approximation is used as a model. The boundary- 
value solution for a circumferential slot operating in the TE^ Q mode 
mounted on a circular cylinder exists and it will be used for comparison. 

The two-dimensional diffraction geometry is shown in Figure l4. It 
was found by comparison with the boundary-value solution that one approxi- 
mation was to replace the edge 1 geometry by a half -plane (n = 2) for the 
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-Diffracted-field (multiple diffractions! 
-Diffracted-field (single-double diffractions] 



270 ° 


180 ° 


Figure 15. Field plot of axial infinite slot on circular conducting 
cylinder with slanted feed (TEM mode). 
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first order diffraction and by a finite wedge (n 0 2) for second and 
higher order diffractions for the 0°-l80° pattern measured in the 
counterclockwise direction. The edge 2 geometry was replaced by a 
finite wedge for all orders of diffraction. For the l80°-360° pattern, 
the approximations of edge 1 are valid for edge 2 and vice-versa. 

The wedge-diffracted field is given by 


V0) = 


1 3 (k|r/2)sin0+acos0cos p] - n/k) f -jkwsinjZfl 


(81) 


_ (1) (h) 

R 1 " R D1 R D1 


(82) 


r 2 . 


(85) 


where R^, R^, ^02^ are gi ven "by Equations 48, 51, 

49, and 52 with 0^ - 90°, Rp^(-0) - 0 and = 0. The creeping 

wave contribution is given by 


E.<0) = E ,e- t l ( “°B +JIc) + E,e‘ t 2 (a - +3k) 


‘B 


(84) 


where 


R A ~ E h 



(85) 
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R B = Ep 



as given by Equation 8l ‘ and 


a = — (ka) 

os 2a 


1/5 1 


9rt 

JlI 


,2/3 


a */6 


( 86 ) 


(87) 


The regions where the appropriate wedge-diffracted and creeping wave 
fields are applicable are identical to those shown in Table 1 for the 
TEM mode. 

The computed results using this approximation are shown on Plate VT 
along with the boundary-value solution. It is noted that the solution 
which includes higher order diffractions gives the best results. A 
noticeable discontinuity is present at the shadow boundary when higher 
order diffractions are neglected. 

Another approximation for the parallel-plate -cylinder geometry would 
be a thin- walled guide (n = 2 for all orders of diffraction) as was 
also true for the TEM mode. Computed results using this model are shown 
in Figure 15* A good agreement is indicated. The discontinuity in the 
} = ±90° directions for the single-double diffractions is again present 
since higher order diffractions are neglected. 

El evation Plane Pattern of Circumferential Slots on Finite Length 
Cylinders Operating in the TEM Mode c 

The elevation plane pattern {j(9)|0 = o3 for a length 

conducting cylinder can be computed using wedge diffraction techniques. 
The diffracted field from the ends of the cylinder and their contribution 
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Diffracted-field (multiple diffractions! 

Diffracted-field (single-double diffractions) 

Boundary-value 



270 ° 


fca = 8.0 


180 ° 


Figure 15* Radiation pattern of circumferential slot on circular 
conducting cylinder using half -planes (TE]_q mode). 
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to the overall radiation pattern can be readily computed. While 
boundary- value solutions for circular conducting cylinders of infinite 
length exist, 3 no solution which takes into account finite length is 
available. To check the validity of the technique, experimental results 
were used for comparison. 

The two-dimensional geometry of a slotted cylinder operating in the 
TEM mode is shown in Figure l6. The diffracted field for the TEM mode 
from edges 1 and 2 is given by 


Ed^o^o) 



( 88 ) 


where 


R D ( 0 o) = e 


0 (kw/2)sin0 


| R l( 0 o ) + M 0 o ) e 


- jkwsinSQ 


(89) 


and Fp(0 o ) and R 2 (9 0 ) are given by Equations 37 and kl. Diffractions 
from the aperture in the 0 Q = ±jt/2 directions mil be diffracted by 
wedges 3 and ^ which in turn mil be diffracted by wedges 1, 2, 3, and 
6. This process continues to higher orders of diffraction. 

To calculate the single-diffracted field by wedge 3, the reciprocity 
principle (see Appendix) mil be applied. This leads to 



\f 2itkr^ 






( 90 ) 
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where 

r d3^ 6 3^ = % (!) [ v b( z 2 i + Q y J1 ^ (91) 

The single-diffracted ray from wedge 3 in the direction of wedge 5 
( e 5 = 180°) w ill he diffracted again and its contribution is 



r c , 0 c- 

y 5 


) = 



( 92 ) 


where 


R^ ) (9 5 ) = R^ ) ( J c)(v B (a,e 5 ,n 5 )} (95) 

Diffractions from wedge 5 in the directions of wedge 3 (©5 =0°) and 
wedge 6 (0^ = 2Y0 o )' will cause additional diffractions. Assuming that 
the length of the cylinder is several wavelengths, the diffractions by 
wedge 6 are negligible. However, the second order diffractions from 
wedge 3 are given by. 


E D3^ r 3^ e 3^ " 


-3 jkr^ +rt/10 




(9*0 


where 

^ ) (0 5 ) = 45 )(9 5 = °)(V a ^ “ W) 


( 95 ) 
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The same procedure used to calculate the fields diffracted by 
wedges 3 and 5 can he applied to calculate the fields diffracted from 
wedges 4- and 6 being initiated by the aperture diffractions in the 


0 Q = -jt/2 direction. 

The first order diffracted fields by wedges 4 and 6 are given by 


(ii -j(K + */3 , . 

%jij. ( r i^®4^ " i %)4 ( e 4) 

y 2itkrij. 

(96) 

where 


-%(- 2 " 2 - 

(97) 

and 


fll -3(kr 6 +*A] 

^ (r ^> “ ^ ^ 

(98) 


where 


4s >(e 6> - ( v B< a ' at - e 6’ n 6>) (99) 


and the second order diffracted field from wedge 4 by 
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^1^ + Jt/9 

^rckri). 




( 100 ) 


■where 


^f(V = *i6 (2jt) ( V B (a ’« + w) 


,(D 


( 101 ) 


The above derivations were all based upon the validity of the 
reciprocity principle as outlined in the Appendix. 

Shifting the coordinates from the individual edges to the center 
of the cylinder and assuming the far-field approximations 




r 0 « r - — cos 6 
° 2 


r^ « r - d cos (7 - 0) 


r, » r + d cos(rt - 7 - 9) ) 


r » r - d cos(it - 7 - 9). 


Ts « r + d cos(7 - 0) 


For phase terms 


( 102 ) 


r « r, » r., « r c » x c r 

o 5 A 5 6 


For amplitude terms (103) 


9 « 0 « 0 , « e « 0 , » 0 
0 3^56 


(104) 


where 
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( 105 ) 


7 = 



(106) 


the diffracted fields, suppressing the e -y \! 2jtkr factor, 

are shown in Table 2 valid in the indicated regions. 

It should be pointed out that additional diffractions by each 
wedge will be present and can very conveniently be accounted for. 
However, for cylinders with several wavelength - radius and length such 
dif fractions have secondary effects and can be neglected. It is 
apparent from Figure l6 that the two-dimensional geometry used for the 
formulation of the elevation plane pattern of a finite length slotted 
•cylinder is identical to the cross section configuration of a slotted 
. finite width ground plane. 

The computed results using the fields from Table 2 are shown on 
Plate VII where they are compared with experimental data.- It is obvious 
that a very good agreement between theory and experiment is indicated. 
The ripples in the -sr/2 < 9 < jr/2 region are present because of the 
diffraction contribution from the edges of the cylinder, and they 
become negligible as the length becomes large. The experimental results 
were obtained from slotted ground plane models since it was deduced . 



TABLE 2. WEDGE-DIFFRACTED FIELDS FOR THE ELEVATION PLANE PATTERN 
OF A FINITE LENGTH CYLINDER IN THE DIFFERENT REGIONS 



rola a ■ roia 
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that their two-dimensional geometry is identical to that of a slotted 
cylinder. 

The beamwidth of the main lobe is decreased as the aperture width 
is increased. The smoothest pattern and the lowest back lobes are 
obtained for the aperture whose width is exactly one wavelength since 
complete cancellation of - the fields diffracted from wedges 1 and 2 in 
the- 0 = ±n/2 directions occurs. The ripples in the -90° < 0 < 90° 
region which are much in evidence in (a) and (b) do not appear in (c) 
and (d) because the fields in the 0 = ±jt/2 directions are very weak 
and any diffractions from the edges are negligible. Additional lobes 
appear as the. aperture width is larger than one wavelength. 

Another formulation for the diffractions from wedges 1 and 2 

37 

would be to use continuous double diffractions. In this case, double- 

— — 

diffracted rays from wedge 2 ^ ( 0 ) are neglected in the 0° < 8 < 90° 

region while double-diffracted rays from wedge 1 Rpi (6) are neglected 
in the -90° < 0 < 0°- region. Computed results using the above formula- 
tion are shown on Plate VIII where they are compared with experimental 
results. Again a very good agreement between theory and experiment 


Is indicated, 



Plate ¥111. Elevation plane patterns of finite length cylinder using 
continuous double diffractions (TEM mode). 




CHAPTER IV 


ELLIPTICAL CYLINDER RADIATION 

Flush-mounted antennas axe commonly used on space vehicles, 
missiles, and- aircraft. However, their radiation properties usually 
cannot he predicted analytically because the geometric shape of the 
body to which they are mounted does not conform to a coordinate system 
where the wave equation is separable. The method used to analyze the 
radiation properties of axial slots on circular conducting cylinders 
can be applied to cylindrical shaped bodies of arbitrary convex cross 
section. One geometry of wide interest is an elliptical cylinder which 
can be used as an approximate model for the wings of an aircraft or the 
fins of a missile. 

The far-field radiation -produced by a slot of arbitrary shape on 
the surface of an elliptical .cylinder of infinite length using modal 
solutions has been carried out by Wait. ^ Computations using the modal 
solution are numerically convenient only for very thin-shaped elliptical 
cylinders, and as the size increases (larger ka and kb) higher order 
terms must also be included for convergence. However, the analysis 
using wedge diffraction and creeping wave theory can be used for com- 
putations of any size and shape of an elliptical cylinder. As the 
physical dimensions increase, the accuracy of the obtained data increases 
because diffraction coefficients (D ^ ) and decay constants ( a^ ) used in 
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the analysis axe derived from asymptotic series of canonical boundary- 
value problems for A 0. 


Equatorial Plane Pattern of Axial Slots Operating in the TEM Mode 

As in the analysis of the circular cylinder, the two-dimensional 
geometry for an axial slot mounted on an elliptical - cylinder and 
operating in the TEM mode is shown in Figure 17. Again, the junctions 
formed by the walls of the parallel-plate waveguide and the planes 
tangent to the surface of the cylinder at the edge points are approxi- 
mated by a pair of infinite wedges of finite included angle 
WA = (2 - n)jt. The wedge-diffracted field is given by Equation 66 and 
is repeated here. 


EjAo^o) = 


-A {k(r 0 -fa/2)sin0] +it/k } 
\J 2irkr 0 


%(0o) +R 2^o) e 


■jkwsin^" 


(107) 


Using the far-field approximations 

r Q « r - T 0 cos 0 cos [3 For phase terms (108) 


r 


o 


r 


For amplitude terras 


(109) 




(HO) 


Equation 107 reduces to 





- ( 0) , _±_ e 0 ^K|p/a/ S in0+T o cos0cosp] -*A) r H (|j) 

^ ' \/2*k \/r l. 


+ E 2 (0)e- 3kwsin?! J 


where 


= \/ x § + yo 

(112) 

_ + W 
2 

(113) 


y 0 = — \j^ a 2 - w 2 (H 1 *-) 

2a 

with R-j_(0) and R 2 (0) as defined previously. 

The creeping wave field is given by Equation <$5 and repeated here. 


E c (s,0) 



B oh( F l> 

W 


-jkCt-L+s-L) 




+ R, 


B oh( p g) = 

D oh l 'A' l/®2 


rtg 

-Jk(t2+S2) “J c£ oh(p) ( 3- £ 


(115) 


From the geometry of Figure 18, one may obtain the far-field approxima- 


tions 




7 6 


8 


2 


it - 


(H7) 


=2 5:3 r - r cos 8-j 


Sg 555 r + r cos 


S ! * s 2 ~ r 


Equation 115 reduces to 


For phase terms 


(118) 


For amplitude terms (119) 


E P (r,0) = V A ^-4^ e" 5kt le" ° 


a oh (p)ds 




n tc 


e-A® -jfcfc 2 "Jo 
+ V-o e e 


a oh (p)ds 
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(120) 


where 


t -e yh) 

A A D oh(y ) 


jkrcosS^ 


( 121 ) 


V B 


= R 


B 


Poh(P2) 

W 


-jkrcosS^ 

e 


( 122 ) 



77 





(125) 


E 


B = 


N/r 




2 

= E i(Q2^[?oh^i] 


(l2t) 


as obtained from Equation 111. 

Expressions for the radius of curvature and differential arc 
length for an ellipse can be obtained by applying the basic geometric 


ds 


\J(dx ) 2 + (dy ) 2 


fl25) 


P 


3/2 

fi + (y* ) 2 J 
y"! 


(126) 


x = t sin 0 


(127) 


y = r cos 0 (128) 


Using Equations 125, 127 > and 128 lead to 


dx = t cos 0 d0 + sin 0 dr 


(129) 
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dy = -T sin 0 d0 + cos 0 dr 


(130) 


-j. 

ds = Qr 2 (d0) 2 + (dr) 2 ] 


y * = } sin_ 

\a / cos 


_b_ (a 2 cos 2 0 + b 2 sin 2 g 
aP cos^0 


Substituting Equations 127 and 128 into the equation for an ellipse in 
rectangular coordinates 


x 2 y 2 
a 2 b 2 


lead to 


r 2 2 j 2 . 2^7 L// 

[a cos 0 + b sin 0J 


dr = 


ab(a 2 - b 2 )sin 0 cos 0 ^ 


r 2 2 . 2 2 

[a cos 0 + b sin i 


(136) 
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Substituting Equations 132 and 133 into 126, and Equations 135 and- 136 
into 131 yield 


P 


ab 


a^cos 2 0 + b^sin 2 0 
a 2 cos 2 0 + b 2 sin 2 0 


3/2 


(137 •) 


as = ab (a 4 cos 2 0 + b^sin^) 1 ^ ^ (138) 

(a 2 cos 2 0 +. b 2 sin 2 0) ^ 

The arc length and complex attenuation factor of the creeping wave on 
the elliptic surface can' then be expressed as 


t =■ 


12 


ds = ab 




L 


#12 (a^co S 2g( + b^sin20) 1//2 
^11 (a 2 cos 2 / + b 2 sin 2 0 ) ^ 


d0 


(139) 



e j #/6 x 




ab d0 



(ihO) 


where 

011 - trapping point angle of the creeping wave 

012 = exit point angle of the creeping wave 
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The points of trapping, Q 1 and Qg, are given by Equations 113 and 
UA. To determine the points of reradiation, P-l and P 2 , of each 
creeping wave on the surface of the elliptical cylinder, the vector 
cross product will be utilized. ^ As shown in Figure 18 for the 
observation angle 0 (assuming far-field observations) 

n x n x = 0 (1^1) 

where 

n = -x cos 0 + y sin 0 = unit vector perpendicular to the (1^2) 

line from the origin of the 
cylinder to the observation 
point 

unit vector normal to the (l^3) 

cylinder surface at the 

point of tangency 

Applying Equation lhl and. using the equation for an ellipse in 
rectangular coordinates 



x 


2 



(ikk) 


leads to 



The wedge-diffracted and creeping 
e -jkr/ f ac tor, are given by 


Reradiation point (i^5) 


Reradiation point P 2 (l 1 ^) 


(1^7) 


( 148 ) 


fields, suppressing the 
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Ep(0) = 


J'( k &/ 2 ) sin0+T o cos0cosffj - 1 1/4^ 


\/2jtk 


jp a + R 2 e-j kwsin ^J (lif-9) 


E c (0) = V A e 


■jkt. 


la oh(p) d£ 


-jtt 2 

+ V B e 


rt 2 

-J <Wp) d£ 


(150) 




ds 



In general, tie total field is given by 


E t (0) = Ed( 0) + E c (0) 


(15D 


( 152 ) 


(153) 


■where the appropriate fields must be considered in each region as shown 
in Table 3* 

Radiation patterns computed using the fields from Table 3 in their 
respective regions are shown on Plates IK, X, and Figure 19- Experimental 
models were constructed and the obtained data were used for comparison 
since boundary-value solutions are numerically convenient only for very 
thin elliptical cylinder. Harrington^ - points out that the equatorial 
radiation pattern (0 = 90 ° plane) for a circular cylinder with finite 
axial slot is identical to that with' an infinite axial slot. It is 
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TABLE 5. WEDGE-DIFFRACTED AND CREEPING WAVE FIELDS FOR THE EQUATORIAL PLANE PATTERN OF AN ELLIFTICAL CYLINDER IN THE DIFFERENT REGIONS 


Regions 
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< 5* <| - 
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< 1 > C 
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Plate X. Comparison of radiation patterns for smaller 
elliptical cylinder (TEM mode). 


size 
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Figure 19. Radiation pattern of a thin elliptical cylinder and a 
finite size ground plane (TEMmode). 



assumed that such a relation holds for an elliptical cylinder, and the 
experimental results were obtained using finite axial slots. The 
experimental models were also of finite length of about 30 wavelengths ; 
diffractions from the edges of the cylinders should then be negligible. 
Dimension construction accuracies were within ±A/25 acd- surface 
irregularities within ±A/50. All measurements were performed at a 
frequency of about 10 GHz. 

The agreement between computed and experimental patterns is good 
particularly in the forward hemisphere. Errors of the order of 2-3 dB 
are obtained in the penumbra region and they indicate that further 
Improvement could be made on the approximations of the parallel-plate - 
cylinder geometry by a wedge and those implicit, in Equations 123 and 
12k* The use of the half -plane wall for certain diffractions processes 
would yield improved results near the penumbra region as has been 
discussed previously for the circular cylinder. In general, the errors 
in the shadow region are of the same order of magnitude, and they are 
quite acceptable for determining possible interference and noise that 
might be introduced from other systems through these back lobes. 

For elliptical cylinders, with .constant major axis but decreasing 
minor axis (thinner shape ellipse), the rate of attenuation of the 
creeping- waves, whose main contribution is in the "shadow" region, is 
larger. However, the "lit" region radiation, whose primary source is 
the wedge-diffracted fields, is essentially unaffected as shown on 
Plate IX (a) -(c). The creeping waves are strongly influenced by the 
radius of curvature which is rapidly changing while the wedge-diffracted 
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fields depend on the included angle of the wedges (WA) which remains 
essentially constant. 

For the small elliptical cylinder (ka = 10.666, kb = 8) of 
Plate X(a) slightly larger variations and asymmetries between experi- 
ment and theory are noted. For ±1.0 percent changes in frequency 
(Plate X(b), (c) ), the variations and asymmetries are even greater and 
it is believed that the physical construction tolerances of the ellip- 
tical cross section are critical for the measurements on this size 
body. For a thin elliptical cylinder the radius of curvature of the 
surface near the major axis is small and any diffractions from that 
region would not be very accurate since the asymptotic series assumption 
of large physical dimensions is not well satisfied. 

In Figure 19, the experimental pattern of a thin elliptical 
cylinder (ka = 40, kb = 8) is compared with a theoretical curve of a 
finite size ground plane obtained using wedge diffraction techniques. 

A good agreement is noted in the 90° < 0 < -90° region. However, in 
the 90° > 0 > -90° region the ripples predicted by the computed curve 
of the ground plane do not appear in the experimental curve for the 
thin elliptical cylinder. The surface of the cylinder is smooth and 
has no sharp edges to contribute diffractions in the "lit" region for 
the ripples to appear. The fact that edge diffracted fields do disturb 
the pattern is an interesting fact to the antenna designer who would 
certainly prefer the smooth pattern. Cylindrical caps on the sides of 
the ground plane would reduce the ripples and smooth the pattern. 

To check the validity of the half -plane (n = 2) approximation for 
the parallel-plate -cylinder geometry instead by a finite wedge (n 0 2), 
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Equations 14-9 and 150 were used with modified n's and the computed 
results were compared with experimental results. It was found by com- 
parison, like for the circular cylinder, that the best approximation was 
to replace the edge 1' geometry by a half -plane (n = 2) for the first 
order diffraction and by a finite wedge (n f 2) for second and higher 
order diffractions for the 0°-l80° pattern measured in' the counter- 
clockwise direction. The edge 2 geometry was replaced by a finite 
wedge for all orders of diffraction. For the l80°-j560° pattern, the 
approximations of edge 1 geometry are valid for edge 2 and vice-versa. 

The computed results using the above approximations are shown on 
Plate XI, where they are compared with the experimental curves. It is 
seen that the above approximations result in better accuracies in the 
penumbra region than the finite wedge approximation but larger varia- 
tions exist in all other regions. As concluded previously, the solution 
used then depends on the region of space of interest. It should be 
noted that computed results for single-double diffractions and multiple 
diffraction shown on Elates IX, X, and XI are almost the same, and it is 
rather difficult to distinguish any variations between the two". 

Equatorial Plane Pattern of Circumferential Slots Operating in the 
TEpo Mode 

As was explained previously, the TE 10 mode is represented by 
two plane TEM mode waves reflecting obliquely back and forth between 
the waveguide walls with the electric field being parallel to the edge 
of the walls. The finite wedge approximation for the parallel-plate- 
cylinder geometry will not satisfy the boundary conditions and the 
half -plane approximation will he used. 
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For a circumferential slot operating in the TE^q mode and mounted 
on a circular cylinder it was found that the "best approximation was to 
replace edge 1 "by a half -plane (n = 2) for the first order diffraction 
and hy a finite wedge (n j- 2) for the second and higher order diffrac- 
tions for the 0°-l80° pattern measured in the counterclockwise direction. 
The edge 2 geometry was replaced by a finite wedge for all orders of 
diffraction. For the l80°-360° pattern, the approximations of edge 1 
geometry are valid for edge 2 and vice-versa. 

The above approximation will be used also for the elliptical 
cylinder. The wedge- diffracted field is given by Equations 81-83 with 
n' s modified as explained above. The creeping wave field is given by 


E C (0) 



(15*0 


where 



eP */ 6 x 


X 


02 


ab d0 


^ [(a^cos 2 ^ + b^sin 2 ^) (a 2 eos 2 0 + b 2 sin 2 0f] 


l7i 


(155) 


The regions where the appropriate wedge-diffracted and creeping wave 
fields are applicable are identical to those shown in Table 3 for the 


TEM mode. 
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The computed results for the TE-^q mode slot on an elliptical 
cylinder are shown on Plate XII. Like for the circular cylinder, the 
single-double diffraction solution has a discontinuity in the 0 = ±90° 
directions because higher order diffractions are neglected. The thin- 
walled guide approximation (n = 2 for all order of diffraction) could 
also be used as was demonstrated for the circular cylinder. Good 
results would also be obtained in this case. 



s sr 




CHAPTER V 


MUTUAL COUPLING BETWEEN SLOTS ON A CYLINDER 

The coupling between two or more waveguide apertures on the 

surface of a cylinder is calculated using wedge diffraction and creeping 

30 

wave theory. Mutual coupling is defined as follows: 


Mutual Coupling = 


Modal current (or voltage) in coupled guide 
Modal current (or voltage) in transmitting guide 


( 156 ) 


It should be. noted that this quantity gives the phase as well as the 
magnitude of the coupling, and may be readily converted to mutual 
impedance or any other desired form for expressing coupling. The 
analysis of coupling will be formulated in the following manner. The 
transmitting guide, with its associated modal current, is replaced by 
an equivalent line source. The modal current induced in the coupled 
guide by a line source is determined. These two steps are then combined 
to obtain the mutual coupling. 

The coupling problem considered in this discussion is described in 
Figure 20. The formulation of the solution to the coupling problem is 
as follows. A unit -amplitude wave with its associated modal current is 
incident in the transmitting guide. An equivalent line source, having 
an omnidirectional pattern and a field matching that of the guide in 
the direction 0^, is substituted for the transmitting guide. A coupled 


9b 
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wave with its associated modal current is induced in the receiving guide. 
The mutual coupling is defined as the ratio of the modal quantity in the 
coupled guide to that of the transmitting guide. 


TEM Mode 

For the TEM mode, the field component perpendicular to the plane 
of Figure 20 is the magnetic field. The incident field is uniform 
across the guide width, and thus the incident power flow for a unit- 
amplitude magnetic field is given by 


P 0 — wZ 0 , 



(157) 


The modal current associated with this mode is 

I 0 = ft • (158) 

The replacement of the guide by a uniform equivalent line source is 
illustrated in Figure 21. The field radiated by the equivalent line 
source is equated to the field radiated by the guide, thus 


-j(kr - it/b) p-jhr 

= H T0 ^ r ^t^ = D T0^t^~^ 


(159) 


where % 0 (^t ) coefficient of the transmitting guide 

in the direction 0^. The modal current 1^ is that of the equivalent 
line source with an antenna impedance Z 0 , and is given by 


-5 


(160) 


I 1 = c T0^t^ e 




A principal criterion for the validity of this analysis is the accuracy 
of substituting the equivalent line source for the guide. 


The modal current induced in the receiving guide by the equivalent 
line source is obtained by application of reciprocity, as illustrated 
by Figure 22. The modal current induced in a line source by a guide 
excited by a modal current 1-^ will be obtained. The power received 
by the line source, being the product of the effective aperture width 
of the line source l/2it, multiplied by the power density, is given by 



- % 2 (Qg)J 



( 161 ) 


where Hrp-^(Q^) and 
points Q-^ and as shown in Figure 22. The modal current received 

by the line source is thus given by 


H T2 (%! ^ axe ^- e fields of the receiving guide at 
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By reciprocity, this • is the modal current induced in the guide by the 
line source excited by a modal current Ip. 

The mutual coupling for the TEM mode is obtained by taking the 
ratio of the modal current in the coupled guide to the modal current 
in the transmitting guide. Using Equations 15 8 , 160, and 162, the mutual 
coupling is given by 



(165 ) 


Computed results of mutual coupling between two slots operating in 
the TEM mode on a circular conducting cylinder are shown in Figure 25 
for different size cylinders and compared with data from a flat ground 
plane. ^ It is seen that the coupling decreases as the relative position 
between the two slots Increases and it would be zero when ■ 
the two slots are diametrica lly opposite since the two creeping "waves 
would cancel in the second slot. It should be noted that the phase of 
the coupling is also available from this analysis. Since theoretical 
determination of coupling by conventional methods is difficult and 
experimental data are not available, no comparison can be made. 



Mutual coupling, 



CHAPTER VI 


CONCLUSIONS 

Wedge diffraction and. creeping wave* theory have been used, to 
analyze the radiation properties of slots on cylindrical bodies'. Using 
this technique, the computed results for circular and elliptical 
cylinders have been favorably compared with existing modal solutions and 
experimental results. Three models representing the parallel-plate- 
cylinder geometry have been introduced and their validity verified. 

One model representing the parallel-plate -cylinder geometry was a 
set of Infinite wedges with finite included angle (n ^ 2). Such an 
approximation served as a good model for TEM mode slots. Although the 
physical structure of the antenna was comprised of flat and curved 
surfaces for ming an edge, a model of two flat surfaces forming an edge 
was a good approximation for computing the diffracted fields in the lit 
region. The values of the fields obtained in the penumbra region may 
not have been very accurate, because of the presence of the wedge, but 
it was demonstrated that good values can be obtained. Such a model can 
be used in approximating other complex structures which may be similar 
to the ones applied in this dissertation. 

Another model representing the parallel-plate -cylinder geometry for 
all orders of diffraction was a set of two flat ground planes (n = 2). 
Even though the physical structure of such a model was quite different 
from the actual antenna, it was demonstrated that it also served as a 


102 



103 


good model for computing the diffracted fields in the lit region. Since 
wedge-diffracted fields are strongly dependent on the tip of the wedge, 
the edge of the flat ground plane was similar to that formed by the 
plate of the guide and the surface of the cylinder. For this case, the 
fields in the penumbra region were more accurate than for the previous 
model because there was no physical structure to interfere with the 
fields. This model can be used for TEM and TEpg modes and will give 
accurate results, as a first order of approximation, for a variety of 
structures which form .a tip even though the physical structure may be 
different. 

The third model used for approximating the edge region of the 
antenna was a half -plane in combination with a finite wedge in the 
following manner: replace the edge 1 geometry by a half -plane (n = 2) 

for the first order diffraction and by a finite wedge (n ^ 2) for second 
and higher order diffractions for the 0°-l80 o pattern measured in the 
counterclockwise direction. The edge 2 geometry is replaced by a finite 
wedge (n / 2) for all orders of diffraction. For the l80°-360° pattern, 
the approximations of edge 1 were valid for edge 2 and vice-versa. It 
was demonstrated that this model also served as a good replacement for 
the parallel-plate -cylinder operating in the TEM and TE^q, modes. It can 
also he used as a model for the wedge-diffracted fields of structures’ 
forming edges which may be similar to the ones for which it was applied 
in this dissertation. 

As a result of the preceding, the antenna designer now has at his 
disposal three models for approximating edges foimed by curved and flat 
surfaces. Each one will give good results as a first order 



of approximation. However, since each one is more accurate in different 
regions, a combination of them will give better results. 

For the TEM mode slots, all three models have been used in the 
computations. The computed results using the half -planes model and 
the half -plane in combination with the finite wedges approximation are 
more accurate for smaller guide widths in the penumbra region but the 
finite wedge approximation gives better results in all other space. 
However, as the guide width is increased, the finite wedge approximation 
gives better results in all regions. 

For the TEpo mode slots, the parallel-plate-cylinder geometry 
cannot be approximated by a set of wedges of finite included angle 
because the tangential E-field boundary conditions are not satisfied. 
Therefore it was necsssary to use either the half -planes model or the 
half -plane in combination with the finite wedges approximation as a 
model. The computations agreed very closely with existing boundary- 
value solutions. 

The elevation plane pattern of a finite length cylinder of an 
arbitrary cross section was computed using wedge diffraction techniques. 
The diffractions from the edges of the cylinder are taken into account 
and their contributions to the overall pattern are noted. As the 
aperture width is increased, the beamwidth of the main lobe is decreased. 
The smoothest pattern and the lowest back lobes are obtained for the 
guide whose width is exactly .one wavelength since complete cancellation 
of the fields diffracted from wedges 1 and 2 along the cylinder 
surface (0 = ±90°) occurs. As the guide width is increased beyond 
one wavelength, additional lobes appear. 
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The maximum radiation of a slotted cylinder can be oriented in any 
desired direction by simply slanting the input guide. Such an orienta- 
tion may be warranted to achieve desired relation characteristics. The 
problem of mutual coupling between slots- on a cylinder has been out- 
lined. Computations have been carried out but cannot be compared since 
boundary- value or experimental data are not available. 

The above technique used to analyze the radiation properties of 
aperture antennas on cylinders presents two problems: approximation 

of the aperture geometry by wedges and the excitation mechanism of the 
creeping waves. In other words, the representation of a complex 
geometry by simpler models whose solutions are known and the initiation 
of the creeping waves seem to be the main obstacles. It was found that 
a Junction formed by a flat plate and a curved surface can be approxi- 
mated by a wedge with infinite sides but finite included angle (n ^ 2), 
by a half -plane (n = 2) or by a half -plane in combination with a finite 
wedge as far as the diffracted fields are concerned. This agrees with 
Soramerf eld’ s conclusion that patterns on precise diffraction photographs 
exhibit almost no dependence on the shape of the diffraction edge and 
the fields diffracted by the wedges in the "lit"' region are chiefly 
dependent on the region near the tip of the wedge. 

All models give good results as has been demonstrated. Also it 
was found that the creeping waves can be initiated by using the wedge- 
diffracted field along the boundary which separates the "lit" and 
"shadow" regions as their initial value. This type of assumption 
preserves continuity of fields as it is necessary. However, this cannot 
very easily be established for any case. For example, creeping wave 
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theory along with image theory on curved surfaces was applied for the 
analysis of the radiation properties of a vertical infinitesimal stub on 
cylindrical bodies. The results obtained were compared with existing 
boundary-value solutions but a very poor agreement was indicated. A 
more practical case of a finite dipole might yield a creeping wave 
solution by segmenting the dipole. Each segment of the dipole would 
launch a creeping wave on the cylinder at the tangent point fixed by 
the height of the segment. There would also be a different "image'' 
element for every radiation angle. This would represent a lengthy 
problem but its solution should be tractable. 

TEM mode aperture field distributions can be realized in practice 
using H-plane sectoral horns of several wavelengths. This type of field 
realization was used for the elevation plane pattern experimental 
measurements . The computed results agreed very closely with the experi- 
mental data which indicates that the assumed aperture field distribution 
is closely approached. However, for circular cylinders, the equatorial 
plane radiation pattern function of an infinite axial TEM mode slot is 
the same as that of finite slot. It is assumed that such a relation 
holds for elliptical cylinders and experimental measurements of TEM mode 
axial slots were carried out using standard rectangular waveguide feeds. 

A comparison between computed and experimental results indicated a very 
good agreement which Justifies its validity. 

The essential feature of this technique is that it is applicable to 
all types of complicated problems some of which are, in general, 
impractical to solve rigorously. The approach is to resolve a complicated 
problem into simpler component problems with approximations, if necessary. 
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such that each of them will have a simple rigorous solution. The known 
solutions of these individual components permit the calculation of the 
overall radiation pattern by the use of superposition. Other analyses, 
such as the point -matching technique have been developed for the solu- 
tion of complicated antennas. However, the point -matching approach 
yields good results for antennas and supporting structures whose 
electrical size is small, resulting in a set of U linear equations 
which are within the capacity of present digital computers to handle 
and the time involved in such computations is great resulting in a high 
cost per data point. The analysis presented in this dissertation based 
on the geometrical theory of diffraction and creeping wave theory yields 
approximate results with- relatively little expenditure of computer time. 
Also, since the field contribution from each individual component is 
separated, this technique is ideal in parametric design problems. 
Individual parameters can be varied and their effect to the overall 
pattern studied. 

Another possible approach to the solution of a complicated problem 
would be an asymptotic series expansion for each problem of interest. 
However, such an approach presents- difficulties. First of all, an 
asymptotic series expansion may not always be very convenient due to 
the complexity of the formulation. Secondly, - even if such an analysis 
is possible, it is applicable only to that individual problem at hand. 

A slight variation in the geometry would require a new; derivation which 
in some cases may not be possible. The technique outlined in the main 
text, although less accurate, provides solutions to general type of 
problems with only slight variations. It should be pointed out that 
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data obtained with the above technique are well within measured 
accuracies. 

The concept of a hybrid solution utilizing wedge diffraction and 
creeping wave theory has been established in solving aperture antenna 
problems on cylindrical bodies that have not been considered previously 
from this point of view or problems that otherwise may not be possible 
to solve in any other way. The first attempt was to apply this tech- 
nique to bodies with known boundary-value solutions for comparison of 
the results. It was then extended to more complicated geometries where 
exact solutions are not possible but where experimental data are avail- 
able. Once a good understanding of the radiation mechanism is 
established, the technique can be extended to other more complex 
geometries to include slots on spheres, prolate spheroids, ogives, 
etc. Also problems involving dipoles oh the surface of such geometric 
shapes should be considered. Experimental data and modern computational 
methods such as the point-matching boundary-value techniques should be 
used to obtain improved values and gain better understanding of the 
radiation mechanisms involved. 

Another possibility would be to combine solutions of the type 
given in this dissertation with the computer techniques. For example, 
the case of the monopole antenna on a cylinder would use point-matching 
methods to find the fields in the vicinity of the monopole and the 
diffraction concepts would be used over the remainder of the cylinder. 
This would greatly extend the size cylinder that could be treated using 
point-matching techniques. 
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DIFFRACTION BY A PERFECTLY CONDUCTING WEDGE 

The diffracted field for a plane wave incident on a wedge of 
included angle (2 - n)it shown in Figure 2k is given by 

Ep(r,t,^ 0 ,n) = V B (r,i|r - t Q ,n) ± V B (r,t + f Q ,n) 


(166) 


where the plus sign applies for the polarization of the electric field 
perpendicular to the edge 


5E 

5n 


= 0 


(167) 


wedge t 


and the minus sign -applies for polarization parallel to the edge 


I wedge j 0 . 

The wedge diffraction function [v B (r # i|r + \]r Q , n )] for a plane wave has 
been determined by Pauli‘S and improved by Hutchins.^ 

The improved form of the diffraction function is defined as 


(168) 


V B (r,t + t Q ,n) = + * 0 ,n) + + t Q ,n) 


( 169 ) 


where 
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eJkra / e^^dr 
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+ Higher order terms 


(170) 


a = 1 + eosQf + t 0 ) - SnidTj (171) 

and N is a positive or negative integer or zero which most nearly 
satisfies the equation 

2m®T - (f + t Q ) = -x for 

2 nnN - (t + t Q ) = +it for 

Equation 170 contains the leading terms plus higher order terms 
which are negligible for large values of hr. For large values of kra 
Equation 170 reduces to the form presented by Pauli and given by 



V B (r,f + t Q ,n) 


a -j(kr+ n/b) 
\f 2 idsr 


I sin (§) 


cos 


St 
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cos 


ft + tc 


n 


(173) 


The diffracted field of Equation 173 is that from which the asymptotic 
diffraction coefficients of the geometrical theory of diffraction are 
obtained.^ This expression is not valid in the shadow boundary because 
a = 0 there. 

The solution for cylindrical wave diffraction at large distances 
from the edge can be determined by the use of the principle of 

c 

reciprocity together with the solution for plane wave diffraction. 
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For plane wave incidence the diffracted field at observation point P 
of cylindrical coordinates (r,t) as shorn in Figure 25(a) is given by 
■V_(r,i,n) of Equation 169. Wow consider the situation in Figure 25 (b) 
in which the wedge is illuminated by a cylindrical wave with its source 
at (x 0 ,g 0 ). Ey reciprocity the diffracted field V A in the direction 
g is equal to the diffracted field V B which is located at the point 
(r = x D ,f = g Q ) with a plane wave incident from the direction = i. 
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(b) 


Figure 25. Illustration of reciprocity. 


BIBLIOGRAPHY 


1. Silver, S. , and Saunders, W. K. "The Radiation from a Transverse 
Rectangular Slot in a Circular Cylinder," Journal of Applied 
Physics, 21 (August, 1950 745-749* 

2. Bailin, L. L. "The Radiation Field Produced by a Slot in a Large 
Circular Cylinder," IRE Transactions, AP-3( July, 1955 ) > 128-137* 

3* Wait, J. R. "Radiation Characteristics of Axial Slots on 
Conducting Cylinder," Wireless Engineer, (December, 1955 )> 

316-322. 

if-. Harrington, R. F. Time Harmonic Fields . New York: McGraw-Hill, 
1961, Pp. 245-254. 

5 . Wait, J. R. Electromagnetic Radiation from Cylindrical Structures . 
Hew York: Pergamon Press, 1959., Pp- 88-104. 

6 . Rudduck, R. C. "Application of Wedge Diffraction To Antenna 

Theory," Report 1691-13, Antenna Laboratory, The Ohio State 
University Research Foundation; prepared under Grant NsG-448, 
National Aeronautics and Space Administration, Office of Grants 
and Research Contracts, Washington, D.C., June 30, 1965* 

7- Ryan, C. E. , and Rudduck, R. C. "Radiation Patterns of Rectangular 
Waveguides," IEEE Transactions on Antennas and Propagation, AP-I 6 
(July, 1968 ), 488-489. 

8 . Ryan, C. E., and Rudduck, R. C. "A Wedge Diffraction Analysis of 
the Radiation Patterns of Parallel-Plate Waveguides," IEEE 
Transactions on Antennas and Propagation, AP-l 6 (July, 1968 ), 

49 O- 49 I. 

9 . Rudduck, R. C., and Yu, J., S. "Higher-Order Diffraction Concept 
Applied to Parallel-Plate Waveguide Patterns," Report I 69 I-I 6 , 
Antenna Laboratory, The Ohio State University Research Foundation; 
prepared under Grant NsG-448, National Aeronautics and Space 
Administration, Washington, D.C., October 15, 1965 . 

Keller, J. B. "Diffraction by an Aperture," Journal of Applied 
Physics, 28 (April, 1957 ) > 426-444. 


114 


10 . 



11 . 


115 


Keller, J. B. "Geometrical Theory of Diffraction," Journal of 
Optical Society of America, 52 (February, 1962 ), 116-130. 

12. Franz, W. , and Depperman, K. "Theorie der Bengang Am Zylinder 
unter Berucksichtigung der Kreichwelle, " Ann. Physik, 10 (June, 

1952), 361-373. 

13* Levy, B. R., and Keller, J. B. "Diffraction by a Smooth Object," 
Communications on Pure and Applied Mathematics, 12(1959), 159-209. 

14. Keller, J. B., and Levy, B. R. "Decay Exponents and Diffraction 
Coefficients for Surface Waves on Surfaces of Won-Constant 
Curvature," IRE Transactions, AP- 3 (Dec ember, 1959 552-561. 

15 . Levy, B. R. "Diffraction by an Elliptic Cylinder, " Journal of 
Mathematics and Mechanics, 9 (i960), 147 - 165 . 

16. Kouyoumjian, R. G. "Asymptotic High-Frequency Methods," Proceedings 
of the IEEE, ( 1965 ), 864-875. 

17 . Sommerfeld, Arnold. Optics . New York: Academic Press, 195^ 

Pp. 245 - 265 . 

18 . Pauli, W. "On Asymptotic Series for Functions in the Theory of 
Diffraction of Light," Physical Review, 54 (Dec ember, 1938), 924-931. 

19. Hutchins, D. L. "Asymptotic Series Describing the Diffraction of a 
Plane Wave by a Two-Dimensional Wedge of Arbitrary Angle," Ph.D. 
Dissertation, The Ohio State University, Electrical Engineering 

D epartment , ( 1967 ) * 

20. Oberhettinger, F. "On Asymptotic Series for Functions Occurring 
in the Theory of Diffraction of Waves by Wedges," J. Math, and 

. Phys., 34(1955), 245-255. 

21. Senior, T. B. A., and Goodrich, R. F. "Scattering by a Sphere," 
Proc. IEE, 111 (May, 1964). 

22. Hong, S. "Asymptotic Theory of Diffraction by Smooth Convex 
Surfaces of Nonconstant Curvature," Technical Report No. 2, The 
University of Michigan, Department of Electrical Engineering 
Radiation Laboratory, Contract AF 04(694) -834, August, 1966 . 

23. Kinber, B. Ye. "Asymptotic Solution of the Problem of Diffraction 
by a Sphere," Telecommunications and Radio Engineering, No. 10, 
October, 1966 . 

24. Ryan, C. E. Jr. "Memorandum on Analysis of Echo Area of Targets 
Using Geometrical Theory of Diffraction and Creeping Wave Theory," 
Report 2430-1, ElectroScience Laboratory, The Ohio State University, 
Department of Electrical Engineering, Columbus, Ohio, May 22, I 967 . 



n6 

25- Peters, L. Jr., and Ryan, C. E., Jr. "The Relation of Creeping 
Wave Phenomena to the Shadow-Zone Geometry," GISAT III 
Symposium, 1967. 

26. Ryan, C. E. "A Geometrical Theory of Diffraction Analysis of the 
Radar Cross Section of a Sectionally Continuous Second-Degree 
Surface of Revolution," Technical Report 2430-4, The Ohio State 
University, ElectroScience Laboratory, Department of Electrical 
Engineering, March, 1968. 

27. Kinber, B. Ye. "Short-Wave Asymptotic Diffraction of Acoustic and 
Electromagnetic Waves from Surfaces of Revolution," Radiotechnika, 

20(1965), 101. 

28. Russo, P. M. , Rudduck, R. C., and Peters, L. Jr. "A Method for 
Computing E-Plane Patterns of Horn Antennas," IEEE Transactions 
on Antennas and Propagation, AP- 13 (Mar eh, 1965), 219-224. 

29. Obha, Y. "On the Radiation Pattern of a Comer Reflector Finite in 
Width," 1 h:h! k Transaction on Antennas and Propagation, AP-ll(March, 
1963), 127-132. 

30. Dybdal, R. B. , Rudduck, R. C., and Tsai, L. L. "Mutual Coupling 
Between TEM and TE-, q Parallel-Plate Waveguide Apertures," IEEE 
Transactions on Antennas and Propagation, AP-14 (September, 1966), 
574-580. 

31. Mikuteit, S. "Mutual Coupling in a Three -Element, Parallel-Plate 
Waveguide Array by Wedge Diffraction and Surface Integration 
Techniques," Report 2485-1, The Ohio State University Electro- 
Science Laboratory, Columbus, Ohio, August 30, 1967* 

32. Rudduck, R. C., and Tsai, L. L. "Aperture Reflection Coefficient 
of TEM and TEqo Mode Parallel-Plate Waveguides," IEEE Transactions 
on Antennas and Propagation, AP-16 (January, 1968), 83-89. 

33. Burnside, W. D. "The Reflection Coefficient of a TEM Mode Symmetric 
Parallel-Plate Waveguide Illuminating a Lossless Dielectric Layer," 
Technical Report 1691-25, The Ohio State University, ElectroScience 
Laboratory, May, 1968. 

34. Jones, J. E., Tsai, L. L., Rudduck, R. C., Swift, C. T. , and 
Burnside, W. D. "The Admittance of a Parallel-Plate Waveguide 
Aperture Illuminating a Metal Sheet, " IE EE Transactions on 
Antennas and Propagation, AP-l6( Sept ember, 1968), 528-535* 

35. Kraus, J. D. Electromagnetics . Hew York: McGraw-Hill, 1953, 
pp. 446-450. 

36. Wu, D. C. F., Rudduck, R. C., and Pelton, E., "Application of a 
Surface Integration Technique to Parallel-Plate Waveguide Radiation 
Analysis," IEEE Transactions .on Antennas and Propagation, AP-17 
(May, 1969 ) • 



117 


37* Wu, D. C. 3?., and Rudduek, R. C., "Slope Diffraction Analysis of 
TEM Parallel-Plate Guide Patterns, " Technical Report 1691-29, The 
Ohio State University, ElectroScience Laboratory, Department of 
Electrical Engineering, Columbus, Ohio. 



